Realization of a classical counterpart of a scalable design for adiabatic quantum 

computation. 
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We implement a classical counterpart of a scalable design for adiabatic quantum computation. 
The key element of this design is a coupler providing controllable coupling between two bistable 
elements (in our case superconducting rings with a single Josephson junction playing the role of a 
classical counterpart of superconducting flux qubits) . The coupler is also a superconducting ring with 
a single Josephson junction that operates in the nonhysteretic mode with a screening parameter of 
about 0.9. The flux-coupling between two bistable rings can be controlled by changing the magnetic 
flux through the coupler. Since the coupling amplitude is proportional to the derivative of the 
coupler's current-flux relation, the coupling can be tuned from ferromagnetic to anti-ferromagnetic. 
In between the coupling can also be switched off. 

PACS numbers: 85.25.Dq 



The magnetic properties of a single-junction interfer- 
ometer (superconducting loop with one Josephson junc- 
tion) depend on its normalized critical current — 
2-kLIc /3>o only. Here L is the loop inductance, Ic the 
critical current of the junction and <E>o is the flux quan- 
tum. If f3 > 1 such an interferometer exhibits a double 
degenerated energy state if an external flux equal to half 
a flux quantum is applied (degeneracy point )i. Close to 
the degeneracy point, the single junction ring is a bistable 
element with two magnetic moments corresponding to 
the superconducting screening currents flowing clockwise 
and counterclockwise. These two states can be described 
in spin formalism by making use of Pauli matrices. In 
other words, a system of magnetically coupled interfer- 
ometers is a realization of a two-dimensional Ising spin 
system. When the coupling between spins is randomly 
distributed the system represents a spin glass. The prob- 
lem of finding the ground state of such a system is a non- 
polynomial one, i.e., the amount of calculation resources 
needed grows exponentially with the number of elements. 
Mathematically this task is equivalent to solving the so- 
called MAXCUT problem^ 

Recently, the implementation of an adiabatic quantum 
algorithm 3 by making use of superconducting flux qubits 
has been proposed^^ Schematically, the implementation 
is the same as for the 2D Ising model described above 
but instead of classical bits quantum bits (qubits) are 
used. It has been shown that such a quantum system 
could solve the problem discussed above in polynomial 
time. In this letter we present the realization of a scalable 
classical design implementing the 2D Ising model with 
a coupling which can be tuned from antiferromagnetic, 
through zero, to ferromagnetic. 

Theoretically several schemes of tunable coupling were 
proposed^ and also some implemented;^ In the design 
presented here we use the approach proposed by Maassen 
van den Brink et al^ Here a coupler, realized as a single- 



junction interferometer with fico = 2ttLcoIc/^o < lj 
provides tunable coupling between two bistable elements. 
The coupling amplitude J can be varied by changing the 
magnetic flux through the coupler. The value of J is 
proportional to the derivative of the coupler's current- 
flux relation 7 s ($ e ), which presents itself as a periodic 
function versus external magnetic flux 4> e . Therefore 
both antiferromagnetic and ferromagnetic couplings can 
be realized in situ. The coupling between interferome- 
ters could also be switched off when dl s /d<& e = 0. We 
have experimentally demonstrated both mentioned types 
of coupling and have been able to switch the coupling off 
in between. 

The scalable design is build up out of basic cells, 
each containing one bistable element directly coupled to 
a readout DC-SQUID (superconducting quantum inter- 
ference device) and two equal couplers coupled by flux 
transformers to the bistable element. In the basic cell 
there are also flux bias lines for each coupler and one for 
both the DC-SQUID and the bistable element. For the 
demonstration of the basic working principle we used a 
structure containing two basic cells of this scalable clas- 
sical design. The samples were fabricated using our stan- 
dard technology with Nb-AlO^-Nb junctions^ The cir- 
cuit diagram and a photograph of this investigated struc- 
ture are shown in Fig[TJ 

The bistable elements (denoted by q, I and q, r, see 
FigOJ are designed as interferometers with f3 q ~ 2.2. 
The coupler is also designed as an interferometer with 
(3co ~ 0.9. Two bistable elements, are coupled via a 
coupler and two flux transformers. The latter are used 
in order to increase the coupling between q, I and q, r 
and separate them in space. This resulted in reduction 
of unwanted crosstalk between the flux bias lines. Each 
bistable element shares a part of its inductance with the 
DC-SQUID used for its readout. The value of this in- 
ductance is optimized in such a way that the biasing of 
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FIG. 1: (a) Photograph of the sample, (b) Circuit diagram. 
Two bistable elements (q,l and q,r) are coupled to each other 
by the coupler. The signals from bistable elements are read 
out by left and right SQUIDs. The fluxes through the inter- 
ferometers could be controlled by the currents (ib.lj ^b,D Ico) 
through the appropriate bias lines (b,l; b,r; and b,CO). All 
junctions are nominally the same: 3.5 x 3.5 /xm 2 . 



the DC-SQUID does not change the state of the bistable 
element. A flux bias line coupled to this shared induc- 
tance can be used for the adjustment of the external flux 
through the DC-SQUID and interferometer as as well as 
for feedback in a flux-locked-loop regime. All Josephson 
junctions are nominally the same (3.5 x 3.5 /im 2 ). The 
desired screening factors for each element are achieved 
only by the loop's geometry. The effective inductances 
of the bistable elements, couplers, and DC-SQUIDs are 
designed taking into account the screening due to the 
coupling with closed superconducting loops. Their de- 
sign values correspond to: L q — 60 pH, Lqo = 25 pH, 
-ksQUiD = 50 pH. 

The mutual inductances between bistable elements, 
DC-SQUIDs, couplers and appropriate bias lines were 
obtained from the measurements of switching current 
histograms as a function of the current through these 
lines. The mutual inductances found from these mea- 
surements are: M bi i- S QuiD,i ~ M b>r]S QuiD, r = 5.4 pH, 

M b,l;q,l ~ Mb,r;q x r = 16 pH and M(, iC 0;Coupler = 6.4 pH. 

For the characterization of the coupler we measured 
the current induced in the right bistable element as a 
function of the current circulating in the left one, i.e. 
the response of the right DC- SQUID to the left bias-line 
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FIG. 2: Demonstration of tunable coupling between two 
bistable elements, (a) The output signal of the left DC- 
SQUID versus the current Ib, T applied to the right bias line. 
The curves correspond to fixed bias current h.co = 524, 526, 
528, 530, 532, 534, 536, 538 and 540; 600, 700 and 800; 850, 
852, 854, 856, 858, 860, 862 and 864 ^A (from top to bot- 
tom) producing dominant magnetic flux in the coupler. The 
readout DC-SQUID operates in flux-locked-loop mode, (b) 
Curves from the upper panel for coupler's bias line current 
h,co =526, 532 and 540 fiA (from top to bottom), the first 
curve is close to the transition from antiferromagnetic to ferro- 
magnetic coupling. Since the magnetic flux through the cou- 
pler weakly depends also on h,r, the gradual crossover from 
antiferromagnetic to ferromagnetic regime is driven by the 
right bias line current. The up, up-down, and down arrows de- 
note antiferromagnetic, "zero", and ferromagnetic couplings, 
respectively. 



current (see Fig[T]) and vice versa. For the measurements 
the DC-SQUID was shunted with an external resistor sol- 
dered onto the sample holder. The DC-SQUID was op- 
erated in flux-locked-loop mode making use of standard 
room temperature electronics^ 2 -. The electronics ampli- 
fies the DC-SQUID voltage and supplies the feedback 
current through the appropriate bias line, keeping the 
total flux in the DC-SQUID constant. For our design 
it means that the output signal of the DC-SQUID elec- 
tronics is proportional to the screening current in the 
interferometer which is roughly proportional to the flux 
coupled to it because q > 1. The measurements were 



carried out in liquid He at 4.2 K. 
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FIG. 3: Output signal of the right DC-SQUID versus the 
current If, t i applied to the left bias line. The current trough 
the coupler are from top to bottom Ib,co = 530, 528, 526, 
524, 522, 520 and 518; 510, 460, 410, 360 and 310; 210, 208, 
206, 204, 202, 200, 198 and 196 ,uA. 

The coupling strength J between the interferometers 
q, I and q, r was tuned by applying a magnetic flux <J> e 
through the coupler using its bias line oc Ib,co (see FigH]). 
The value of J increases with the coupler's screening fac- 
tor /3co- On the other hand, /3co < 1 is required in order 
to obtain the change of the coupling sign. Thus, for opti- 
mal results we need a sample with /?co ^$ 1- Taking into 
account the coupler inductance Leo — 25 pH, we fab- 
ricated samples with the junction oxidation process for 
j c = 100 A /cm 2 giving a critical current for the junction 
Ic ~ 12 fiA and flco ~ 0.9 close to unity. 

The results are shown in Figj2]were the left DC-SQUID 
is used as a measuring device and in Fig [3] for measure- 
ments with the right DC-SQUID. The slope of the curves 
is due to the direct parasitic coupling between the right 
(left) bias line and the left (right) bistable element. The 
periodical modulation of these curves is due to the jumps 
of screening current of the right (left) interferometer. 
Note that the period of this modulation is the same as 
the period measured with the DC-SQUIDs directly cou- 



pled to the interferometer confirming that we are looking 
at the effect of this interferometer. 

In Figs|2]and[3]the curves with upward and downward 
kinks correspond to antiferromagnetic and ferromagnetic 
couplings, respectively. The rounding is due to an av- 
eraging of the thermally activated jumps between two 
stable branches of the interferometer. The kink height 
is proportional to the derivative of the couplers' current 
phase relation: dl/d$ e = 2nlc cos </)/[$ (l + flco cos</>)] 
with <f> = 27r$ e /$o — /3cosin</> the phasedrop over the 
coupler junction. So the maximal kink heights in the 
ferro- and antiferromagnetic cases are proportional to 
/3co/(l ~ Pco) and /3co/(l + Pco), respectively. The 
observed ratio of the maximal kink height for these two 
cases 12:1 corresponds to fico = 0.85 which is in reason- 
able agreement with the expected value. Also, the fact 
that the ferromagnetic coupling region is smaller than 
the anti-ferromagnetic one qualitatively agrees with this 
picture. However, the experimentally observed width of 
this region is only 0.04 $o while this simple model would 
predict 0.23 $o for /3co = 0.85. This could be due to a 
non-sinusoidal current phase relation of the coupler junc- 
tion. Further experiments would be necessary to clear up 
this disagreement. In between the ferro- and the antifer- 
romagnetic regions the coupling vanishes as indicated by 
the updown arrow. 

In conclusion, we demonstrate tunable coupling be- 
tween two single junction interferometers in classi- 
cal mode. We show both ferromagnetic and anti- 
ferromagnetic type of coupling, including zero coupling 
between these two regimes. This system enables to study 
relaxation and annealing in spin glass. 
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